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1. Introduction 

In the last few years several experimental groups have observed the crossover [1] from the 
weakly paired Bardeen-Cooper-Schrieffer (BCS) state to the Bose-Einstein condensation 
(BEC) of molecular dimers with ultra-cold two-hyperfme-component Fermi vapors of 40 K 
atoms [2-4] and 6 Li atoms [5,6]. The unitarity limit of the Fermi superfluid was attained 
by manipulating an external background magnetic field near a Feshbach resonance which 
allows an experimental realization of infinitely large value of the s-wave scattering length 
a [7]. 

The most interesting feature of the BCS-BEC crossover, noted in experiments [8] 
on a BCS superfluid as well as demonstrated in theoretical model calculations [9] is that, 
due to a dominance of the Pauli repulsion among fermionic atoms over the interatomic 
attraction, the Fermi superfluid remains essentially repulsive in the unitarity limit, 
i.e. the gas does not collapse and its properties are quite regular. Elaborated Monte 
Carlo calculations have confirmed this effect [10]. A similar conclusion follows from an 
examination of the compressibility of a Fermi gas [11]. The Fermi system should exhibit 
universal behavior in the unitarity limit which should limit the maximum attractive 
force to a finite value as a — ► ±oo [8]. This phenomenon has greatly enhanced the 
interest in the theoretical study of a Fermi gas in the unitarity limit [12]. 

Close to the critical temperature a Fermi superfluid can be studied by using the 
Ginzburg-Landau theory with a complex order parameter [13-15]. Recently a nonlinear 
Schrddinger (NLS) equation for the complex order parameter has been proposed at 
zero-temperature to study the BCS-BEC crossover [16-19]. This equation has a 
simple nonlinear term in the weak-coupling BCS limit (a negative and small), at 
unitarity (a = ±oo), and the BEC regime (a positive and small) where it becomes 
the Gross-Pitaevskii equation for Bose-condensed molecules. However, in the full BCS- 
BEC crossover the 3D NLS equation has a complicated nonlinear term [7,10,17,20]. 
In addition, this 3D NLS equation must be solved numerically [16-18]. Hence an 
effective one-dimensional (ID) and two-dimensional (2D) reduction of this equation 
under convenient trapping conditions of cigar and disk-shaped traps is well appreciated. 
The present paper addresses this important question of ID and 2D reduction of the 
original 3D NLS equation at unitarity. Through numerical studies we also investigate 
the validity of the 3D- ID and 3D-2D reduction of the 3D NLS equation. 

In Sec. [2] we consider a confined 3D Fermi superfluid at unitarity and discuss the 
3D NLS equation [18,21]. We derive it from an energy functional which takes into 
account the bulk properties of the system and also the inhomogeneities in the density 
profile due to the external potential. We show that the gradient term of our energy 
functional is practically the same of the one recently deduced [22] with an epsilon 
expansion of energy density in 4 — e dimensions [22,23]. In addition, we find that 
the total energy calculated from this equation for a 3D system of spin 1/2 fermions, 
trapped in a spherically-symmetric harmonic potential, is in good agreement [24] with 
those obtained from accurate Monte Carlo calculations [25,26]. In Sec. [3] we introduce 
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a generic axially-symmetric harmonic potential which models the confining trap of the 
superfluid fermionic system. In Sec. H] we consider a cigar-shaped Fermi superfluid. 
Using a variational ansatz we minimize the 3D energy functional and derive an effective 
ID nonpolynomial Schrddinger (ID NPS) equation for this system. Simple analytic 
forms of the ID NPS model are derived for small and large number of atoms. The 
variation of the nonlinear term from small to large number of atoms is illustrated. In Sec. 
[5] we consider a disk-shaped Fermi superfluid. Using a variational ansatz we minimize 
the 3D energy functional and derive an effective 2D nonpolynomial Schrddinger (2D 
NPS) equation for this system. Simple analytic forms of the 2D NPS model are derived 
for small and large number of atoms. In Sec. [6] we present a numerical study of the 
different model equations in one and two dimensions for cigar- and disk-shaped systems 
and compare the results with the full 3D NLS equation. The results of the approximate 
models are found to be in good agreement with exact calculations. Finally, in Sec. [7] 
we present concluding remarks. 

2. Fermi superfluid at unitarity 

We consider a dilute Fermi gas of N atoms with two equally populated spin components 
and attractive inter-atomic strength. At zero temperature the gas is fully superfluid and 
the superfluid density coincides with the total density. A description of the superfluid 
state can be obtained by using a complex order parameter \I/(r) conveniently normalized 
to the total number of superfluid pairs N p [13-15], i.e. 



The local number density of superfluid pairs is defined as n p (r) = |\[>(r)| 2 = n(r)/2 where 
n(r) is the total density (the density of atoms). Notice that \& is not the condensate 
wave function, because the modulus square of it gives the number of particles, not of 
condensate particles [18,21]. At unitarity (a = ±oo) one finds the condensate density of 
pairs is n ~ 0.7 n p , where n p is the total density of pairs [27,28]. In the BEC side (a > 0) 
of the BCS-BEC crossover n reaches n p quite rapidly by reducing the positive scattering 
length a, while in the BCS side (a < 0) of the crossover n decreases exponentially to 
zero by reducing the absolute value of the negative scattering length [27,28]. However, 
the phase of ^(r) is exactly the phase of the condensate wave function and its gradient 
gives the superfluid velocity [18,21]. 

Under an external potential U(r) acting on individual atoms, the properties of the 
Fermi superfluid in the BCS-BEC crossover can be described, in the spirit of the density 
functional theory [29-31], by the energy functional 



where U p (r) = 2U(r) is the external potential acting on a pair and U(r) the external 
potential acting on a single atom [16-18]. Here S[n p ] is bulk energy density as a function 




(2.1) 




(2.2) 
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of density of pairs, i.e. the energy density of the uniform system in the BCS-BEC 
crossover [16-18]. 

The last two terms in ( 12.21) correspond to the local density approximation (LDA), 
equivalent to hydrostatics, and the first term involving the gradient corresponds to a 
correction to LDA [22,31]. Including only the LDA term in energy density leads to 
the Thomas-Fermi (TF) approximation. In the gradient term m p = 2m is the mass 
of a pair with m the mass of one atom. This term takes into account corrections to 
the kinetic energy due to spatial variations in the density of the system. For normal 
fermions various authors have proposed different gradient terms [32-36]. For superfluid 
fermions we are using the familiar von Weizsacker [32] term in the case of pairs. 

We notice that the present gradient term is exactly the one that emerges [37] from 
the Bogoliubov-de Gennes equations in the BEC regime of the BCS-BEC crossover. In 
this case the energy functional (I2.2p is nothing but the Gross-Pitaevskii energy functional 
of the Bose-condensed molecules [37,38]. 



In addition, we stress that the gradient term [h 2 /(2m p )] V-v/n p (r) of (12.21) 
is very close to that recently obtained in [22] for the same gradient term in the 
case of a Fermi gas at unitarity, in a rigorous calculation using an e expansion of 
energy density around d = 4 — e spatial dimensions [22, 23]. Their final result 
for this term 0.032/i 2 [Vn(r)] 2 /[mn(r)], quoted in their (53), can be rewritten as 

r -i 2 

0.512(/j, 2 /m p ) Vyn p (r) , if we recall m p = 2m and n(r) = 2|\I/(r)| 2 = 2n p (r) — 
very close to the present gradient term in (12.21) . 

In the BCS-BEC crossover, the bulk equation of state of energy density £{n) as a 
function of density of atoms n of a dilute superfluid depends on the s-wave scattering 
length a of the inter-atomic potential [16-18]. In the unitarity limit, when a — > ±oo, the 
bulk energy density is independent of a and, for simple dimensional reasons [11,39-41], 
must be of the form 

S(n) = \x~n^ = 2 2 ^ X -nf , (2.3) 
5 m 5 m 1 

where \ is a universal coefficient. Thus, in the unitarity limit the bulk chemical potential 

is proportional to that of a non-interacting Fermi gas. For fermions, results of fixed-node 

Monte Carlo calculation give % = (3vr 2 ) 2/3 {/2 with £ = 0.44 [10]. 

Taking into account (12. 3p , the energy functional of the Fermi superfluid reads 

E = f J-^|V*(r)| 2 + 2U{v)\^{v)\ 2 + 2 2 / 3 - X -|*(r)| 10 / 3 ) d 3 r . (2.4) 
J 1 4m 5 m J 

By minimizing (12.41) with constraint (12. ip we obtain the following nonlinear Schrddinger 
equation 



ft 2 ?H 2 
_ZL V 2 + 2U(r) + 2 2 / 3 X — |#(r)| 4 / 3 
4m m 



(r) = 2/i ^(r) , (2.5) 



where the chemical potential /io is fixed by normalization. This is the zero-temperature 
3D NLS equation of superfluid Fermi gas at unitarity. 
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Figure 1. Ground-state energy E (in units of fku) versus N from a solution of the 3D 
NLS equation (solid line). The results of Green-function Monte Carlo (GFMC) [26] and 
fixed-node diffusion Monte Carlo (FNDMC) [25] are shown for a comparison (symbols) 
The local density approximation (LDA), i.e. the Thomas-Fermi model, is also exhibited 
(dashed line). 



In the case of a large number of fermions, in (12.51) the nonlinear term vastly 
dominates over the kinetic energy term containing the Laplacian V 2 , and the numerical 
solution of this equation is thus insensitive to the Laplacian term. In phenomenological 
treatment of Fermi superfluid, this term is often neglected and the resulting model is 
called the TF model (LDA approach), which, however, leads to a nonanalytic solution in 
space variable. The inclusion of the Laplacian term leads to a solution of (12.51) analytic 
in space variable, although lying close to the solution of the TF model for a large number 
of particles. The Laplacian term, however, plays an important role in the case of a small 
number of particles, as we shall show studying the dimensional crossover of the system. 

Now we compare the energies of our density functional (12.41) with those of Monte 
Carlo calculations [25, 26] in the case of a harmonically trapped system with small 
number N of spin 1/2 fermions at unitarity. We solve (12.51) numerically, using the 
Crank-Nicholson method detailed in Sec. El and calculate the total energy E of the 
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system using ( 12.4ft for different N in a spherically-symmetric harmonic trap 

U(r) = ^mu 2 r 2 (2.6) 

with uj the trap frequency. In figure [1] we plot the total energy E in units of hu vs. 
N. For comparison, we also plot the TF energy. The energy of the TF approximation 
(LDA approach) is analytically known to be E(N) = (3N) 4 / 3 y/^/A [42]. From the data 
of Figure [T] we find that the average percentage deviation of the present result from the 
FNDMC data is 3.30%, whereas the same of the LDA result from the FNDMC data is 
9.47%. Very recently it has been shown [43] that the best fit to the fixed-node diffusion 
Monte-Carlo data at unitarity [25] is obtained, in the case of even number of particles, 
by using A = 0.18 in the gradient term Xh 2 \ V\I/| 2 /m of the energy functional (12 Ah . 

The usefulness of the present model equation (12.51) over the commonly-used TF 
approximation (LDA approach) is well appreciated. In fact, the TF approximation is 
valid only for very large N values when the gradient term is negligible compared to 
the bulk chemical potential. The present density-functional model for Fermi superfluid 
at unitarity has been extended over the full BCS-unitarity crossover and the extended 
model has been found to yield [24] energy of a Fermi superfluid under spherical harmonic 
confinement in good agreement with Monte Carlo data [44] in the crossover domain. 



3. External confinement: axially-symmetric harmonic potential 

If we take r = (p, z) so that p and z are the radial and axial variables and if we consider 
the axially-symmetric harmonic trap 

U(v) = ^mco 2 ± (\ 2 2 p 2 + \ 2 1 z 2 ), (3.1) 

where A2^j_ and XiU± are the trap frequencies in radial and axial directions with Ai and 
A2 convenient trap parameters. In the following we shall consider for simplicity a real 
wave function ^(r). Equation (12. 5ft can be reduced to the following dimensionless form 
by scaling z = za±, p = pa±, = tyyj N/2/a 3 / 2 , and p = pohtxi± 



A\dp 2 dz 2 J Apdp 

+ 2N 2 / 3 x^ 4/3 } ^(p, z) = 2p ij(p, z) , (3.2) 



where a± = ^h/mu± is the characteristic harmonic oscillator length in the radial 
direction. Equation ( 13.21) satisfies normalization 2ir J °° dp J^° oo dzp^ 2 (p, z) = 1. 

For a cigar-shaped superfluid it is convenient to define a linear density through 

POO 

f 2 (z) = 2n dpp^ 2 (p,z). (3.3) 
Jo 

Similarly, for a disk-shaped superfluid it is convenient to consider a radial density 
through 

POO 

2 (p)= / dz^ 2 (p,z). (3.4) 
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4. Cigar-shaped Fermi superfluid: 3D-1D crossover 

4-1. Harmonic Confinement 

Now let us suppose that the external trapping potential U(r) is given by a harmonic 
confinement of frequency u± in the cylindrical radial direction p and by a generic 
potential V(z) in the cylindrical axial direction z: 

U{t) = ^mu 2 ± p 2 + V{z) . (4.1) 
We introduce the variational field 

into the fermionic energy functional (12.41) and integrate over x and y coordinates. 
After neglecting the space derivatives of a (z) (adiabatic approximation) we obtain the 
following effective energy functional 



Ei 



1 4m 



df(z) 
dz 



2 



+ 



2V(z) + 



h 2 



Ama 2 (z) 



a 2 (z)] ~f(~z) + 2^^- h - } -^-\ dz (4 3) 

° [Z) \ 1 [Z) + 25tt 2 /3 m o*lHz) J ' 1 J 



which depends on two fields: the transverse width a(z) and the axial wave function 
/>")■ 

Note that the variational approach we are using here has been successfully applied 
in the dimensional reduction from 3D to ID of the 3D Gross-Pitaevskii equation, which 
describes trapped Bose-Einstein condensates. In that case the variational approach ends 
up with the ID nonpolynomial Schrddinger equation (ID NPSE) for the axial wave 
function [45]. The ID NPSE has been extended to investigate the Tonks- Girardeau 
regime [46,47], the two-component BEC [48], transverse spatial modulations [49] and 
also axial vorticity [50], 

Minimizing E\ with respect to f(z) one finds 

^ + 2 ^ ) + 4^M + m ^ 2(5) 

f(z) = Tp x f(z) . (4.4) 



,o2/3 6X h 2 f±l\z 

57r 2 / 3 m a 4 / 3 (z 



This is a ID Schrddinger equation and p\ is fixed by the normalization 

f 2 (z) d~z=^ = N p . (4.5) 
Instead, minimizing E\ with respect to a one gets 

0I V( 2 ~) = i + 2^£^jni)« 2,3 (i> ■ (4.6) 

We call (14.41) . equipped with (14.61) . the ID nonpolynomial Schrddinger (ID NPS) 
equation. 
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The ID NPS equation can be conveniently written in dimensionless form by scaling 
za±, a(z) = a(z)a±, f(z) = ^N/2f(z)/y/a ± , and pt,\ = fku±(fii + 1/2) as follows: 



+N 2/3 6 X f 4 '\z 



f(z) = 2^f{z) , (4.7) 



57r 2/3 a 4/3( z 

a\z) = \ + N^^fr\ z)(T ^{ z ), (4.8) 

with normalization f 2 (z) dz = 1 . Here we have included a constant term [=1/2, 
corresponding to the energy of the system in the transverse trap whose effect has been 
integrated out in (14.71) and (jHEJ) ] in the scaled chemical potential, so that in the N = 
limit the chemical potential coincides with that of the axial harmonic trap. In the N = 
limit, from (14.81) we find that o~ 2 (z) = 1/2 and the 1 / (4cr 2 ) + a 2 terms cancel the constant 
term in (14.71) . In deriving (14.71) we assumed an harmonic confinement in the z direction: 
V(z) = m\\uj 2 _z 2 jl with Ai the anisotropic parameter. 

By using the ID NPS equation (14.71) with (14.81) we can study the dimensional 
crossover from 3D to ID of the superfluid Fermi gas at unitarity. In general (14. 8 j) must 
be solved numerically for a(z) in terms of f(z) and the result when substituted in 
(14 .7j) leads to the principal result of this section. No closed-form analytic expression 
for this equation can be given as in the bosonic case [45], except under special limiting 
conditions. 

The first interesting limit of the formalism is obtained for small number of atoms 
when the nonlinear term in (14.TD is small, (which corresponds to the weak-coupling 
limit in the bosonic case,) so that the last term in ( 14. 8ft can be neglected and the 
transverse width is z independent and is given by a(z) = a = l/y/2 under the condition 
N P f 2 (z) < 125 v / 2tt/(128(3x) 3/2 ) = 0.273293... (obviously satisfied for a small number 
of fermions). The cigar-shaped system is then quasi-lD and governed by 

1 d 2 6iV p 2/3 x(2/) 4 / 3 1 

\~dz~ 2 + V + 5^ J f{z) = 2/il f{z) • ( 9) 

In the opposite extreme, for a large number of fermions, N p f 2 (z) ^> 
125\/27r/(128(3x) 3/2 ), the cigar-shaped system is effectively 3D. Under this condition 
(USD can also be solved for a(z) to yield a(z) = [4\/iV P /(^)/v / 7r] 2/5 (3x/25) 3/10 . 
Consequently, the 1/(4<t 2 (2:)) term in ( 14. 7p can be neglected and the remaining terms 
combined to yield 

Id 2 . 2 2 7 (6^f^i 
4^ + A ^ + 5 (5^)1/5 J / W = 2 ^ • 

(4.10) 

The power of the nonlinear term has changed from 7/3 to 9/5 as we pass from 
quasi-lD regime governed by (14 ,9p to the effectively 3D regime governed by O4.10p . 
In the quasi-lD regime the nonlinear power 7/3 is the same as that in the original 
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Figure 2. (a) The s vs. n\ dependence as calculated from a numerical solution of 
(|4.13p and (|4.14|) . (b) Square of the frequency £1 2 = 2 + s in a cigar-shaped trap 
vs. density n\ as calculated from (|4.15|) . (c) Sound velocity c\ vs. density n% of 
a uniform gas from a numerical solution of (|4.17p and (|4.18p and the asymptotic 
Ci = 0.766846ny 5 result valid for large n\. All results refer to the unitarity limit 
£ = 0.44. 

three-dimensional equation (12 .51) . whereas in the 3D regime it has acquired a different 
power. Although the quasi- ID equation (14. 9 p has been used in different studies on a 
Fermi superfiuid [51] and on a degenerate Fermi gas [52], (I4.10p for Fermi superfiuid is 
new. We shall see that (14. 10p is already valid for a moderate number of Fermi atoms 
(N > 100) and is of interest in phenomeno logical applications. 

In the 3D regime it is a good approximation to neglect the kinetic energy term in 
(14. 10p and the following analytic expression for density is obtained in the so called TF 
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approximation 

N P f(z) = j^f^i - W 2 ©(2^ - A?. 2 ), 

« - A 2 , 2 ) 5 / 2 6(2 / i 1 - A 2 , 2 ) (4.11) 

where 0(x) is the Heaviside step function. As we are in the 3D regime it is interesting 
to compare this result with the following TF approximation made on the full three- 
dimensional equation (12.51) after integrating over the transverse variables 

N P f{z) = ^f-^i - \\z^Q{2^ - \\z% 

« - \\z 2 f'M^i - A?* 2 ). (4.12) 

The two TF results have the same functional dependence on the variables as well as 
very similar numerical coefficients, in spite of (14.101) and (12.51) having different powers 
of density in the nonlinear terms. The quasi-lD equation (14. 9ft has the same power of 
density in the nonlinear term as (12. 5p . Nevertheless, a TF approximation made in (14.91) 
will generate a density with an entirely different dependence on z. 

It is now interesting to study the bulk chemical potential /ii(ni) implicit in (14.71) 
and (14.81) given by 

Mm) = T^T-T + o\z) - 1 + 2 2 / 3 -^^ (4.13) 



4a 2 (z) w 5vr 2 / 3 a 4 /3( z ) 



4 '"^ 3 + 22/3 2ili' ! ' /3(l)ff2/3(l) ' (4 - 14) 

where density n\{z) = N p f 2 (z). Many physical observables are determined by the 
bulk chemical potential so obtained. For example, assuming a power-law dependence 
yt,i ~ n\ for the bulk chemical potential on density n\ (polytropic equation of state), 
the frequency of the lowest axial compressional mode Vt\is given by [46] 

2 



it 



2 + s. (4.15) 

Here we introduce an effective polytropic index s as the logarithmic derivative of the 
bulk chemical potential /ii, that is 

s = —t— ■ (4.16) 
Hi oni 

From (14.131) and (14.141) one finds that in the ID regime s = 2/3 and Q = Oi/(Aiu;_i_) = 
a/8/3, while in the 3D regime s = 2/5 and f2 = Jli/(Aio;j_) = a/12/5. Note that 
the result a/12/5 is exactly the same one obtains setting 7 = 2/3 in the formula 
(37 + 2)/(7 + 1) derived by Cozzini and Stringari [53] for the axial breathing mode 
of 3D cigar-shaped superfluids with a 3D bulk chemical potential that scales as nj. 

In figure [2] (a) we show the power s of the dependence fii ~ as a function of rii 
obtained by numerically solving nonlinear (I4.13P and (I4.14p . In figured (b) we plot the 
square of collective frequency Q 2 in the dimensional crossover as a function of the axial 
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density n\ and As the number of atoms iV increases, the frequency of axial compressional 
mode slightly decreases with a decrease of the polytropic power in \i ~ n\ relation. 

4-2. Uniform Density 

Now we consider the axially uniform case, where V(z) = 0. In this case all the 
variables attain a constant value independent of z and we remove the z dependence 
on all variables. Setting n\ = N p f 2 , from (14. 7\\ and ( 14.81) we find 

2/3 

(4.17) 



1 1 2 1 2> X 



2ni 



ixo 2 



1 , 12x r 



a 



4 25 



2n 1 cr 



7T 



(4.18) 



These equations can be used to derive the axial sound velocity c\ of the fermionic system, 
that is obtained with the formula [14] 



c, = Jn^ . (4,9, 



dp 
dn 

In the quasi-lD regime (for small ni), from (14.171) and (14.181) we obtain a 2 = 1/2, 
pi = 3x(4ni/7r) 2 / 3 /5 and ci = ^2x/5(Ani/n) 1 ^ 3 . In the effective 3D regime 
(for large n 2 ), from (14171) and (14~T8D we obtain a 2 = (12x/25) 3/5 (2rii/7r) 2 / 5 , /ii = 
(7/5)(3x) 3/5 (n 1 / 7 r) 2 / 5 /20 1 / 5 and Cl = v /l4725(3x) 3/10 (ni/7r) 2 / 10 /20 1 / 10 = 0.766846n} /5 
for ^ = 0.44. In figure [5] (c) we plot the sound velocity function of the axial 

density n\ as calculated from a full numerical solution of (14.171) and (14.181) as well as the 

1/5 

asymptotic result C\ = 0.766846ri 1 valid for the effective 3D regime for large n\. The 
two results are indistinguishable except near the origin. This shows that the effective 3D 
description (14.101) is very good except for very small atom number. The sound velocity 
increases with matter density as it should. 

5. Disk-shaped Fermi superfluid: 3D-2D crossover 

Let us suppose that the external trapping potential U (r) is given by a generic potential 
W(p) in the cylindrical radial direction p and by a harmonic confinement of frequency 
u z in the cylindrical axial direction z: 

U(r) = W(p) + l -mu) 2 z ~z 2 . (5.1) 

We introduce the variational field 

1 ( z 2 



^ ] =i^^)^\:ww) ]m (5 - 2) 



into the fermionic energy functional (12.41) and integrate over the z coordinate. After 
neglecting the space derivatives of fj(p) we obtain the following effective energy 
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functional 



2vr 



2W(p) 



h 2 



lmrf(p) 



-mujjfj 2 (p) 



010/3 (ft 



2 (p)+2 2 / 3 



2X 



7T 



1/3 



3/2 



p dp 



(5.3) 



m ?7 2 / 3 (p) / 

which depends on two fields: the axial width fj(p) and the transverse wave function 

Also in this case we observe that the variational approach we are using has been 
successfully applied in the dimensional reduction from 3D to 2D of the 3D Gross- 
Pitaevskii equation. The resulting effective equation has been called 2D nonpolynomial 
Schrddinger equation (2D NPSE) [45]. 

Minimizing E 2 with respect to <fr(p) one finds 



4m 



V 2 + W(p) + 



ft 2 

— — — + -mu z r] p 
imr] 2 (p) 2 



+ 2 2/3 fl 2X ^ 2 4/3 (p) 

V 5 7T 1 / 3 m r/ 2 / 3 (p) 



4>{p) = 2p 2 0(p) 



(5.4) 



This equation is a two-dimensional Schrddinger equation and p 2 is fixed by the 
normalization 



2tt 



<p 2 (p) pdp = ^- . 



Instead, minimizing E 2 with respect to fj(p) one gets 

^ 3/2 

where a 



a. 



V(p) 



I , o2/3 4 X 

4 + 3^/3 



4> 4/3 m 4/3 ( P ) 



(5.5) 



(5.6) 



^J^JJjmo^) is the characteristic harmonic length in the axial direction. We 
call ( 15.41) . together with ( 15.61) . the 2D nonpolynomial Schrddinger (2D NPS) equation. 

The 2D NPS equation can be conveniently written in dimensionless form by scaling 
p = pa z , fj{p) = rj(p)a z , 0(p) = y/N p (f)(p)/a z , and p 2 = hw z {p 2 + 1/4) as follows: 



- 4V 2 + p \ 2 + 



+iv2/ 3 /3^x t!M 
V 5^1/3^2/3^) 



r/ 2 (p) 
<Kp) 



2p 2 0(p) 



^(p) 



4 + 3vrV3 ^ 



S) 



3/2 



4/3 (p)^ /3 (p) 



(5.7) 



(5.? 



with normalization 27r J °° (j) 2 (p)pdp = 1. In deriving (15. 7p and (15. 8p we assumed an 
harmonic confinement in the radial p direction: W(p) = muj 2 \\p 2 / "1 . Again in defining 
the reduced chemical potential p 2 we have removed the zero-point energy corresponding 
to the energy of the axial trap, so that in the N = limit, (15.71) and (15. 8p coincides with 
the corresponding linear harmonic oscillator problem. 
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By using ( 15.7ft with (15.81) we can study the dimensional crossover from 3D to 2D 
of the superfluid Fermi gas at unitarity. First (I5.8P is to be solved numerically for rj(p) 
in terms of 0(p) and the result when substituted in ( 15.71) gives the desired result for 
studying a crossover from 3D to 2D. Closed- form analytic result for these equations is 
possible only under limiting conditions. 

The first interesting limit of the formulation is obtained for a small number of 
atoms when the nonlinear term in (15. 7\i is small, so that the last term in (I5.8P can 
be neglected. Under this condition N p( p 2 < 25 v / 3^(5/3) 1 / 4 /(192x 3/2 ) = 0.148656... 



(obviously satisfied for a small number of fermions) the longitudinal width rj(p) = y 1/2 
is independent of p. The disk-shaped system is then quasi-2D and described by 



T + p2xl + Vf ^" 04/3 l m = 2/i20(p) 



(5.9) 



In the opposite extreme, for a large number of fermions, N p <p 2 ^> 
25v / 37r(5/3) 1/ ' 4 /(192x 3 / 2 ), the disk-shaped system is effectively 3D. Under this condition 
(EE]) can be solved for rj(p) to yield r](p) = 2A p 1/ V /8 3 3/1 V /2 (p)/(5 9/16 vr 1 / 8 ). Substi- 
tuting this result in (15. 7p and neglecting the l/[8?7 2 (p)] term we get 



- T + f% + 57r i/ 451 /B ' \ M = ^(p) • (5-io) 

The power of the nonlinear term has changed from 7/3 to 2 as we pass from quasi-2D 
regime governed by ( 15.91) to the 3D regime governed by ( I5.10p . In the quasi-2D regime 
the nonlinear power 7/3 is the same as in the original three-dimensional equation (12. 5p . 
whereas in the 3D regime it has acquired a different power. Of the models ( 15.91) and 
(I5.10p . (15. 9p was previously considered by others [54], whereas (15.101) is new. However, 
we shall see that (15.101) should have wide phenomenological application for disk-shaped 
superfluid as this form is already effective for N > 100, producing better approximation 
than CTI . 

In the 3D regime if we neglect the kinetic energy term in ( 15.101) and the following 
analytic expression for density is obtained in the TF approximation after a neglect of 
the kinetic energy term 

0.201862 / O , ON 2 „ / n,n, 

« (2p 2 - p 2 X 2 ) &(2p 2 - p 2 X 2 ). (5.11) 

As we are in the 3D regime it is interesting to compare this result with the following 
TF approximation made on the full three-dimensional equation (I2.5P after integrating 
over the longitudinal variable 

AW 2 (P) = (2A* - P 2 A 2 ) 2 Q(2p 2 - p 2 A 2 ), 

« (2/i 2 - P 2 A 2 ) 2 ®{2p 2 - p 2 \ 2 ). (5.12) 
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Figure 3. (a) The s vs. 712 dependence as calculated from a numerical solution of 
(|5.13p — (|5.15|) . (b) Sound velocity c 2 vs. density n 2 of a uniform gas from a numerical 
solution of (|5.16p and (|5.17p and the asymptotic result C2 - 
All results refer to the unitarity limit £ = 0.44. 



0.986212ri2 / for large n 2 . 



The two TF results have the same functional dependence on the variables as well as 
very similar numerical coefficients, in spite of ( I5.10I) and ( 12. 5ft having different powers 
of density in the nonlinear terms. The quasi-2D equation (15. 9 p has the same power of 
density in the nonlinear term as ( 12. 5ft . Nevertheless, a TF approximation made in ( 15.91) 
will generate a density with an entirely different dependence on p. 

We now consider the bulk chemical potential implicit in (I5.7P and (15. 8p 



^2(^2) 

v 4 (p) = 



1 , l„2 (o) _ I , 02/3 /3 2x n 2 /3 (p) 
lrp{p) 2 ,KP) 2 V 5ttV3 t? 2 / 3 (p) 



1 
1 



2 2/3 



4x 
3ttV3 



3/2 



(5.13) 



(5.14) 



where density n 2 (p) = N p <p 2 (p). Again the bulk chemical potential can be considered 
to possess a power- law dependence on density: p, ~ n|, where the numerical coefficient 



s can be extracted from a numerical solution of nonlinear ( 15.131) and ( 15.141) using the 
relation 

n 2 dp 2 



P2 dn 2 ' 



(5.15) 



Effective Nonlinear Schrddinger Equations at Unitarity 



15 



The coefficient s is of interest in the study of physical observables of interest. In the 
quasi-2D regime s = 2/3 whereas in the 3D regime s — 1/2. However, in the quasi- 
2D to 3D crossover the /i(n 2 ) dependence is to be calculated numerically using (15.131) 
and (I5.14p and then the polytropic indexs can be calculated as a function of n 2 for the 
dimensional crossover. In figure [3] (a) we plot the polytropic power s vs. density. 

5.1. Uniform Density 

Now we consider the radially uniform case, where W(p) = 0. In this case all the 
variables attain a constant value independent of p and we remove the p dependence on 
all variables. Setting n 2 = N p (p 2 , from (15.71) and (15. 8p we find 



These equations can be used to derive the radial sound velocity c 2 of the fermionic 
system, that is obtained with the formula [14] 



In the quasi-2D regime (for small n 2 ), from (I5.16P and (15.171) we obtain rj 2 = 1/2, 



/i 2 = 2%A/3/5n 2 /7T 1 / 3 and c 2 = 2Wx/vT5n<j/ 3 /7r 1 / 6 . In the effective 3D regime 



(for large n 2 ), from (l4"T7j) and (OS1) we obtain n 2 = 4( v / 3x) 3/4 r4 / 7(5vr 1/4 5 1/8 ), /i 2 = 
(6/5)(27/5) 1 / 8 x 3/4 n 2 /2 /7r 1/4 and c 2 = v ^75(27/5) 1 / 16 x 3/8 n 2 /2 /vr 1/8 = 0.986212r4 /2 for 
£ = 0.44. In figure [3] (b) we plot the sound velocity c 2 vs. the axial density n 2 as 
calculated from a full numerical solution of (15.161) and (15.171) as well as the asymptotic 

1 /2 

result c 2 = 0.986212^2 for large n 2 in the effective 3D description given by (I5.10p . The 
asymptotic result is indistinguishable from the exact result except near very small n 2 . 
This shows that the effective 3D description of the system is very good. The sound 
velocity increases with matter density as it should. 

6. Numerical Result 

Next we study the effectiveness of the dimensional reduction of the 3D GL equation (12.51) 
to ID and 2D forms with a variation of the number of fermions for cigar- and disk-shaped 
configurations. We numerically solve the full 3D equation as well as various ID and 2D 
reduced equations by discretizing them by the semi-implicit Crank- Nicholson algorithm 
with imaginary time propagation [55-57]. For numerical convenience we transform 
the chemical potential term pip in the nonlinear equations to a time-dependent term 
idip/dt. The space and time steps used in discretization were typically 0.05 and 0.001 
respectively. 




(5.16) 



(5.17) 




(5.18) 
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Figure 4. Normalized density f 2 {z) f 2 (z)dz = 1) along the axial z direction for 

(a) 2, (b) 10, and (c) 100 Fermi atoms from a solution of 3D equation (|3.2|) denoted 3D, 
the full one-dimensional equations (|4.7|) and (|4.8p denoted NPS, the quasi-lD equation 
(|4.9p denoted Ql, and the effectively 3D equation <|4.10[) denoted E3 in the unitarity 
limit using parameters £ = 0.44, Ai = 0.1, and A2 = 1. 

6.1. 3D- ID crossover 

For a cigar-shaped superfTuid we solve four sets of equations in the unitarity limit: (a) 
the 3D equation (13.21) . (b) the complete reduced ID equations (14.71) and (14. 8p . (c) the 
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approximate quasi-lD equation ( 14. 9ft . and (d) the approximate effectively 3D equation 
(14. 10p . We use the parameters £ = 0.44 [10], and a highly cigar-shaped trap with 
Ai = 0.1 and A2 = 1. In figure [4] we illustrate the results of our calculation by plotting 
the linear density profile f 2 {z) vs. z of the four sets of calculations for fermion number 
N = 2, 10 and 100. 

From figureHJwe find that the ID approximate calculations are good approximations 
to the solution of the 3D equation (13.21) . However, some features of the different 
approximate models are worth commenting. All approximate ID models lead to a 
density smaller than that obtained from the 3D equation (13. 2p . The density obtained 
from (14. 71) and (14.81) provide the best approximation to the exact density for all N. For 
small N, the quasi-lD model (14 .9p provides a better approximation to the exact result 
than the effectively 3D model ( 14. 10[) . The opposite happens for large values of N. For an 
intermediate value of N, the quasi- ID model ( 14.91) and the effectively 3D model (14.101) 
could produce similar results. 

6.2. 3D-2D crossover 

For a disk-shaped superfluid we again solve for sets of equations: (a) the 3D equation 
( 13. 2L (b) the complete reduced 2D equations ( 15. 7ft and ( 15. BJ) , (c) the quasi-2D equation 
(15.9p . and (d) the effectively 3D equation 1 15. 10ft . We use the parameters £ = 0.44 in the 
unitarity limit and Ai = 1, and A 2 = 0.1 for a disk-shaped trap. In figure [3] we present 
the results of our calculation by plotting the radial density profile (p 2 (p) vs. p of the 
four sets of calculations for fermion numbers N = 2, 10, and 1000. 

From figure [5] we find that the three 2D approximate models could be good 
approximations to the solution of the 3D equation 113.21) . Again, as in the cigar-shaped 
superfluid, the full 3D model 113.21) produces the largest density profile with the complete 
2D equations 1 15.71) and 115.81) providing the best approximation to it for all values of N. 
For small N, the quasi-2D model 115.91) produces better approximation to the exact 
result than the effectively 3D model 115. 10p . The opposite happens for large N. For an 
intermediate N these two latter approximations could produce similar results. 

An interesting result from our calculations with cigar and disk-shaped superfluid 
is that, for N > 100, the quasi-lD model 114.91) and the quasi-2D model ( 15.91) are 
poorer approximations than the effectively 3D models (14. lUp and (I5.10p . respectively. 
For experimental purpose N = 100 represent a small number atoms. Hence for 
phenomenological applications the effectively 3D models (I4.10p and (15.101) with nonlinear 
terms with power 9/5 and 2 should be used. Note that these powers are different from 
the power 7/3 in the original 3D equation (13 .2[) . The quasi-lD model (14.91) and quasi-2D 
model (15.91) for cigar- and disk-shaped superfluids with nonlinear terms of power 7/3 
effective for a small number of fermions are useful for academic interest. 




(c) 6 12 18 24 



P 

Figure 5. Normalized density </> 2 (p) (2tt J °° pcf) 2 {p)dp = 1) along the radial p direction 
for (a) 2, (b) 10, and (c) 1000 Fermi atoms from a solution of 3D equation (|3.2|) denoted 
3D, the full two-dimensional equations (|5.7p and (|5.8|) denoted NPS, the quasi-2D 
equation (|5.9|) denoted Q2, and the effectively 3D equation (|5.10|) denoted E3 in the 
unitarity limit using parameters £ = 0.44, Ai = 1, and A2 = 0.1. 
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7. Conclusion 



We have suggested a time-independent Schrddinger equation for a Fermi superfluid at 
unitarity [ f)3.2p ] by minimizing its energy functional. This equation can also be derived 
as an Euler-Lagrange equation of an appropriate Lagrangian density. In a cigar-shaped 
superfluid assuming a Gaussian form for the order parameter, and integrating over the 
transverse variables we have derived an effective nonlinear nonpolynomial ID equation 
for the Fermi superfluid at unitarity [ (14 .7ft and (14.81) ]. This complex equation is simplified 
in the limit of small and large atom numbers when it reduces to a nonlinear equation 
with power-law nonlinearity. The equation for small atom number N has the same 
nonlinear structure as the original 3D equation and is called quasi-lD model [ (14.91) ]. 
whereas the equation for large N has a distinct nonlinearity and is called effective-3D 
model [( 14. 101) ] . For phenomenological application the effective-3D model seems quite 
attractive. 

In a disk-shaped superfluid assuming a Gaussian form for the order parameter, and 
integrating over the axial variable we also derived an effective nonlinear nonpolynomial 
2D equation for the Fermi superfluid at unitarity [( 15. 7ft and (" 15.81) ] . This complex 
equation is simplified in the limit of small and large atom numbers when it reduces to 
a nonlinear equation with power-law nonlinearity. The equation for small atom number 
iV has the same nonlinear structure as the original 3D equation and is called quasi- 
2D model [ (15.9j) ] . whereas the equation for large N has a distinct nonlinearity and is 
called effective-3D model [ (15.101) ]. The quasi-2D model has the same nonlinearity as the 
original 3D equation, whereas the effective-3D model produces a different nonlinearity. 
Again the effective-3D model is attractive for phenomenological application producing 
very good results. All the above models have been studied by a numerical solution of 
the model equations. 
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